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Abstract 

An exactly solvable reaction-diffusion model consisting of first-class 
particles in the presence of a single second-class particle is introduced on 
a one-dimensional lattice with periodic boundary condition. The number 
of first-class particles can be changed due to creation and annihilation 
reactions. It is shown that the system undergoes a discontinuous phase 
transition in contrast to the case where the density of the second-class 
particles is finite and the phase transition is continuous. 

One of the most important characteristics of non-equilibrium driven systems is 
that their steady state consist of current of particles or energy. One-dimensional 
reaction-diffusion models are examples of such systems which have attracted 
much attention during last decade p] [2]. Phase transition and shock forma- 
tion in these systems are some of their interesting collecting behaviors. These 
systems have also many applications in different fields of physics and biology. 
During recent years different models of this type have been studied widely and 
interesting results have been obtained. The Asymmetric Simple Exclusion Pro- 
cess (ASEP) is a well known example. In this exactly solvable model, which is 
defined on an open discrete lattice, particle are injected from the left bound- 
ary and extracted from the right boundary while hopping on the lattice to the 
left and to the right randomly. This model has been shown to exhibit non- 
trivial steady-state phenomena such as phase transitions and shock formation 
[3 . In order to study the steady state properties of these shocks different mod- 
els have been proposed. It should be noted that the ASEP is not the only 
one-dimensional out-of-equilibrium system which exhibits shocks. It has been 
shown that there are three families of two-states models in which a factorized 
shock measure is invariant under the time evolution if some constraint on the 
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microscopic reaction rates are fulfilled [3]. In [5] the authors have shown that 
the same phenomenon might also be observed in three-states systems. 
In this paper we study an exactly solvable three-states model with non-conserving 
dynamics on a discrete lattice with a ring geometry. Our model belongs to the 
class of non-conserving driven-diffusive systems where attachment and detach- 
ment of particles are allowed . The study of such models which are variants of 
the ASEP is motivated by the biological transport processes in living systems 
[51 [7] and the denaturation transition in DNA [51 d] . Since our model is based 
on a newly introduced model in [lOj . we will first briefly review the main con- 
cepts and results of this model. In [lOj a non-equilibrium three-species system 
is introduced on a lattice with periodic boundary condition consisting of the 
following reaction processes 

A% — > %A with rate a 
%B — ► B% with rate (5 

AB — > BA with rate 1 (1) 
A% — ► 00 with rate A 
00 — ► A$ with rate A'. 

As can be seen the number of A particles (first-class particles) is not conserved. 
In contrast, the number of B particles (second-class particles) is conserved since 
they only diffuse. It is assumed that in a system with at least one empty site one 
has finite number of second-class particles with the density ps in the presence 
of the first-class particles with fluctuating density. It has been shown that in 
this case a continuous phase transition takes place if the order parameter of the 
system is taken to be the density of empty sites in the system pe- By taking 
a = (3 = 1 and defining u> := A it turns out that pe is zero for u> < uj c while 
it changes linearly as pe = — Pb for > u> c in which ui c = l PB flB ■ The 
current of the second-class is always constant while the particle current of the 
first-class particles is given by different expressions in each phase. For the case 
a 7^ 1 and (3=1 the transition point is obtained to be lo c = pB 1 ^ a ~ ■ The 
density of the empty sites is zero below the transition point while it is given by 
Pe = yuj ~ i A-'if-a a bove this point. For ps ^ the transition is still continu- 
ous. 

In present paper we assume that there exists only a single second-class particle in 
the system which means their density goes to zero in the thermodynamic limit. 
Second-class or tagged particles are usually introduced to study the dynamical 
properties of the shocks in one-dimensional driven-diffusive systems; however, 
one of our major motivations for studying such limiting case is to investigate 
its effects on the critical behavior of the system and compare it with the pre- 
vious case in which the density of the second-class particles is non-zero in the 
thermodynamic limit. As we will see considering this limiting case changes the 
nature of phase transition from a continuous into a discontinuous one. As far as 
we know such observation had not been reported before. Apart from the vast 
applicability of such models in different fields of science (as mentioned above), 
classification of one-dimensional driven-diffusive models which are exactly solv- 
able using the Matrix Product Formalism (MPF) has been of great interests for 
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people in this field (for a recent review see [H]). As we will see the model is 
still exactly solvable using the MPF even in the limiting case ps — > 0. 
In the following we define u := j, and apply the MPF [3] to find the partition 
function of the system. According to the MPF the stationary probability dis- 
tribution function of any configuration C of the system of length L + 1 with a 
single second-class particle at the site L + 1 is given by 



(2) 



in which Xj = E if the site i is empty otherwise = A. The normalization 
factor Z in the denominator of © will be called the partition function of the 
system. By applying the standard MPF the quadratic algebra of the model is 
obtained to be [TO] 

AB = A + B 
AE = ±E 
EB = |e 
E 2 = ^E. 



(3) 



By defining E = ^|V)(W| in which (W\V) = 1 one finds from © 

AB = A + B 

A\V) = 
(W\B=1(W\. 



(4) 



This quadratic algebra has an infinite-dimensional representation given by the 
following matrices and vectors 
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(5) 



in which a 2 = a ^~ 1 . Since the stationary state of the system without vacancies 
is trivial, we consider the partition function of the system with at least one empty 
site which is defined by 



Z = TV [(A + E) L B] - Tr[A L B]. 



(6) 
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Using (JUnj) and after some straightforward calculations we find the following 
exact expression for the partition function of the system 



0(1 +u — a) a 1 — a a 0(1 + w — a)(l — a) 

For a < 1 there is no phase transition. Assuming a > 1 one simply finds the 
following expressions for the partition function of the system in the large L limit 

(8) 



/3(l+w-a)(l 



for a; < a — 1. 




At the transition point the partition function of the system grows like 0(L). 
Taking the density of the empty sites on the lattice given by 

uj d 

p E = lim — — InZ (9) 

L^oa L OUJ 

as the order parameter of the system, we find using |SJ that 

for uj > a — 1 

PE = { (10) 



At the transition point the density of the empty sites is obtained to be pe — 
2(i +u} ) ■ We should note that density of the empty sites for uj < a — 1 drops to 
zero as pe oc O(j^). Discontinuous changes of the density of the empty sites 
Pe in the thermodynamic limit indicates a first-order phase transition in the 
system. As we mentioned earlier, in the case where the number of the second- 
class particles on the lattice is finite the density of the empty sites p e changed 
continuously over the transition point [lOj . 

In order to study the nature of the first-order phase transition one can apply 
the Yang-Lee theory. Recently it has been shown that the classical Yang-Lee 
theory can be applied to the out-of-equilibrium systems to study their phase 
transitions (for a review see 12]). We have calculated the line of the Yang-Lee 
zeros for our model in the complex-w plan and found that they lie on a circle of 
radius a. The center of this circle is at (—1, 0) and intersects the real-w axis at 
Re(to) = a — 1 at an angle § which again implies a first-order phase transition 
at the transition point. The density of the zeros has also been found to be a 
constant all over the circle. 

It is also interesting to calculate the density profile of the first class particles 
on the ring, as seen by the second-class particle, using the MPF. For a system 
with at least one empty site it is given by 

p A (i) = ]-(Tr[(A + E) t A(A + E) L - l - 1 B]~Tr[A L B]) < i < L - 1. (11) 
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It turns out that (fTTj) can be calculated exactly using ([5]) and here are the results 
in the large L limit 



Pa{i) 



i 

l + UI 



„(a+/9-l) 
a(/3+a>) 



1 - ^e" 



-J L " f±1 



for w > a — 1 



for £j < a — 1 



for u> = a — 1 



(12) 



in which the correlation length is given by £ = |ln(i^)| _1 . For w > a — 1 
the lattice is filled by first-class particles of density except just in front of 
the second-class particle where it increases exponentially to 1. In this phase the 
density of empty sites is . For ui < a — 1 the density of first-class particles 
increases exponentially from — to 1 in the bulk of the lattice. The density of 
empty sites in this phase is nearly zero in the thermodynamic limit. As can be 
seen at the transition point the density profile of the particles is linear. This 
is a sign for a shock however since the number of first-class particles is not 
a conserved quantity the shock position fluctuates and therefore the resulting 
profile is linear. This phenomenon has also been observed in the ASEP with 
open boundaries on the first-order phase transition line where the injection and 
extraction rates become equal and smaller than one-half. The sock picture 
will be more clear by calculating the connected two-point function of first-class 
particles. Straightforward calculations result in the following exact expression 
which is valid for i < j and large system length 

(PA(i)PA(j))c ■■= (p A {i)PA{j)) ~ (pA(i))(PA{j)) 

(13) 

= -(rb - (PA(i)))V - (paU)}). 

We have also calculated the current of the first-class particles J a in the steady 
state. In the large L limit we have found that the current of the first-class 
particles does not depend on j3 and is given by 



Ja 



for lu > a — 1 



for lu < a — 1 



for to = a — 1. 



(14) 



V. 2(l+a;) 

On the other hand, the mean speed of the second-class particle defined as 



1 



V = -^{f3Tr[{A + E) L ~ 1 EB] + Tr[(A + E) L ~ L AB] 



L-l 



Tr\A L B] 



(15) 



can also be calculated exactly. It turns out that V is given by the following 
exact expression in the thermodynamic limit 



V = 



(l+u)(J3+u>) 



for lu > a 
for lu < a 



(16) 
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Figure 1: The mean speed of the second-class particle V as a function of u> for 
three values of (3. The length of the system is L = 100 and we have chosen 
a = 6. 



In Figure Q] we have plotted V as a function of oj for three different values of 
(3 on a lattice of length L = 100. For oj < a — 1 the speed of the second-class 
particle is equal to one and does not depend on oj while for to < a — 1, above 
the transition point, for (3 < 1 ((3 > 1) the speed of the second-class particle is 
a decreasing (increasing) function of oj. For (3=1 the speed of the second-class 
particle is always equal to unity. 

The model studied in this paper consists of a single second-class particle in the 
presence of first-class particles with fluctuating density because of creation and 
annihilation of them. Comparing out results with those obtained in [lOj one 
should note that the order of the phase transition is changed from two to one 
when the density of the second-class particles goes to zero. The physical ex- 
planation for such a macroscopic change can be as follows: for the two cases 
Pb = and ps ^ and below the critical point there are empty sites on the 
ring; however, their density goes to zero in the thermodynamic limit. Never- 
theless, it is less probable to find configurations of type A% in the case ps / 
(in comparison to the case pb = 0) because in this case part of the system is 
occupied by the second-class particles. Therefore as we increase oj above the 
critical point it is more probable to create many empty sites in the case pb = 
than the case pb ^ because as explained above we have more configurations 
of type A$ in this case which result in the configuration 00. This means that as 
we increase oj above the critical point we expect many empty sites to be created 
at once in the case ps = in contrast to the case pb ^ where they are being 
created slowly. 

It is also interesting to compare our results with those obtained in [T3] where the 
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same model as (JT|) has been considered except it does not contain the creation 
and annihilation of the first-class particles. It has been shown that the model 
in this case has two phases: a condensate phase and a fluid phase. Obviously 
in our model by fixing the number of the first-class particles the phase in which 
Pe = will be destroyed; however, a condensate phase emerges in which the 
density profile of the first-class particles is no longer linear but an step-function 
and this is in quite agreement with the results in [13] . 

The MPF enables us to solve some of the one-dimensional driven-diffusive sys- 
tems exactly; however, by looking at these models we realize that they might 
have similar quadratic algebras. This means that a quadratic algebra can de- 
scribe different models with different physical properties and critical behaviors. 
Now the question is that whether or not other three-states models defined on 
a ring geometry can be described by ([3]). In fact we have found that there is a 
family of such models which are exactly solvable and share the same quadratic 
algebra [H]. So far two members of this family are introduced and studied in 
[TU] and [H]. 
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